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ABSTRACT 

The theorem of least work was used to investigate 
the stress pattern and bowing, i.e. uniform curvature, 
of an initially straight thin tube subjected to a steady 
temperature distribution independent of the longitudinal 
coordinate of the tube. Simple equations were developed 
for stress, strain, and bowing. They show that if the 
thermal strain is expressed as a Fourier series, the tan- 
gential stress is small and is a function of the first 
harmonic only, the axial stress is a function only of 
the second and higher harmonics, and the bowing is a 


function only of the first harmonic. 
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Nomenclature. 


A 
a,b 


£@ 


To 


T(6) 


u(@) 


Wega 


n 


v(@) 


cross section area of tube wall 


coefficients used to specify axial strain 
(See Equation 4) 


local thermal strain, presumed a function of 6 
only, not of r or Zz. 


coefficients in Fourier expansion of €@) 
Young's modulus of elasticity 


an integer denoting the nth term of a 
Fourier series 


radial coordinate, measured positive outward 
from mid-thickness of tube 


tube radius to mid=-thickness of wall 
tube wall thickness 
local temperature 


base temperature from which thermal strain is 
measured 


strain energy per unit length of tube 


local radial displacement of mid-thickness 
curve of tube 


= coefficients in Fourier expansion of u(@) 


local tangential displacement of mid-thickness 
curve of tube 


= coefficients of Fourier expansion of v(@) 


= tube coordinates (see Figure 1.) 


coefficient of linear thermal expansion 


a 


t 
ei 
elastic strain components in tangential and 


longitudinal directions, respectively 








oa 5,” 


Vy x © 


< s Hh 


Ww 


z 


}+ oe 


4% 


hn 


total strain (elastic plus thermal) in the 
tangential and longitudinal directions 
respectively 


2 14 

fot age (0-1) 
angular coordinate measured clockwise from 
the positive z-axis, (see Figure 1). 
angular position of the plane of bowing 
curvature of bowed tube 
Poisson's ratio 

Ao 


Le ipnk 


2 
[+ ee (n™ |) 


direct stress in the radial, tangential, and 
longitudinal directions, respectively 


shear stress in the rG@-plane 


rva-EY¥ ) 
ffi + vai-¥)} 
Bava 


So CE 


{/- yp? 


/- FE , 
yt Tash ) 


dw 


ag Where w is any representative variable 


6 





2. Introduction 

The problem of stress distribution in a tube with 
temperature variation independent of the axial coordinate 
has been solved under the following two sets of circum- 
stances: (1) an exact solution(2f*in the sense of the 
theory of elasticity for steady-state temperature config- 
uration, and (2) a very elementary solution[5] in the 
framework of elementary strength of materials. Solution 
(1) is of limited usefulness, not only because of its 
limitation to the steady-state configuration but also 
because of its mathematical complexity. Solution (2) was 
developed to investigate technically important cases of 
the bowing of pipes partially filled with a cryogenic 
fluid, such as liquid oxygen. This solution assumes 
that only longitudinal strains and stresses in the pipe 
are an important cause of bowing. This in effect is 
taking Poisson's ratio equal to zero. Presuming that 
Poisson's ratio is zero, while actually it is in the 
neighborhood of 0.3 for most piping materials, led to a 
simple and useful solution for bowing radius and axial 
stress but said nothing about tangential stress. The 
degree of approximation involved was quite unknown; there- 
fore, it was the purpose of the present investigation to 
obtain a more accurate solution than that of Flieder, 


Loria, and Gnith [5] without involving the mathematical com- 


$$ 
Numbers in brackets refer to works listed in the 
bibliography, page 45. 





Then strain energy is differentiated partially with 
respect to each such coefficient, each such partial deri- 
vative is equated to zero, and finally the resulting set 
of equations is solved for the corresponding values of 
the coefficients. 

This procedure, which was used for the present investi- 
gation, gives the displacement which minimizes the strain 
energy (subject to the implied constraints of the assumed 
displacement) and thus corresponds to equilibrium. 

That the theorem of least work can be applied to the 
problem at hand, where temperature changes are involved, 
is not obvious. Argyris[(il and Boley and Weiner [2| have 
discussed, in general terms, the energy theorems and thermal 
stresses where temperature changes and non-linear elasti- 
city apply. Even so it seems appropriate to establish, 
by using the following physical considerations, that the 
theory of least work can be applied to the present problem. 

Consider the structure cut into a number of elemental 
volumes such that the temperature change is uniform in 
each element. A change in temperature causes thermal 
strains but not stresses; however, the elements no longer 
fit together continuously. In order to fit the elements 
together into a continuous structure having an assumed dis- 


placement pattern from the original shape, stresses must 





aa 





be applied which cause elastic strains in accordance with 
the theory of elasticity. Since the stresses are applied 
holding temperature constant, the strain energy can be 
calculated using the equation previously cited, Jaodr 3 
The values of the elastic strains are determined by the 
difference between the total strain, determined from 
geometric considerations of the assumed displacement pat- 
tern, and the thermal strains resulting from known tempera- 
ture changes. The theory of least work now applies and 
states that the strain energy is a minimum when the assumed 


displacement pattern is that corresponding to equilibrium. 





3. Analysis of Stress and Strain. 

In order to proceed with the analysis, certain condi- 
tions must be met and assumptions made. The tube must be 
constructed of an isotropic elastic material and the 
temperature distribution must be a known function of the 
angular coordinate, O , and independent of time and the 
longitudinal coordinate, z. This is the resulting tempera- 
ture distribution when a tube is partially filled with a 
cold liquid. Also no external forces or moments may be 
applied to the tube. The necessary assumptions are as 
follows: 

1. The tube is "thin" and thus % and 7, are small 
and can safely be neglected in the analysis.” 

2. The temperature is constant across the tube wall 
thickness. 

3. Young's modulus of elasticity and Poisson's ratio 
are constant. This is areasonable assumption as shown by 
values for 302 stainless steel in Figuresl and 2. 
Stainless steel is typical of most metals. 


4. Plane normal cross sections remain plane. 


*By neglecting ¢,and tre,equilibrium cannot be satisfied 
and any results are necessarily approximate. A similar 
assumption has been made by others under similar circumstances 
and the solutions obtained were acceptable when compared 
with the exact solutions according to the theory of elasti- 
city. See DenHartog[4] pages 221 to 223 and page 240. The 
latter portion of Appendix III shows the relative magnitudes 
which may be expected of 6, and 7rg. 
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5. In accordance with Saint Venant’s principle the 
stress pattern exists uniformly along the tube except at 


short distances from the ends. 





“400 -300 -200 -100 O 160 200 300 400 


Figure 1. Modulus of Elasticity of AISI 302 
Stainless Steel 





~400  -300  -200 -100. Oo j00. 200 300 400 
TEMPERATURE, °F 


Figure 2. Prorepon'SsPRatic for Aisi 302 
Stainless Steel 


From T. F. Durham, R. M. McClintock, and R. P. 
Reed, "Cryogenic Materials Data Handbook”, 

~~ Bureau of Standards, PB Report 171809, 
1961. 
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A. Symmetric Case. 
Temporarily the temperature distribution will be 
limited to one symmetric about the y-axis. The coordinate 


system used is shown in Figure 3. 





Figure 3.. Typical Section of Tube 
When there is a change in temperature of the tube, a 
point P located on the mid-thickness curve will, in gen- 
eral, be displaced to a new location P. The radial dis- 
placement, u, and the tangential displacement, v, are 


positive in the directions indicated in Figure 4. 





Figure 4. Displacement of Mid-thickness Curve 





The tangential strain at any point can be expressed as the 
strain of the mid-thickness curve plus the strain due to 
bending of the tube wall, which in turn can be expressed 
in terms of the radial and tangential displacements as 
follows €,* tt a as3. : (1) 
For derivation see Appendix I. 

Since the temperature distribution is symmetric about 
9= 0, the radial displacement will be an even function of 


@ and the tangential displacement will be an odd function 


of O@ . Expressed as Fourier series they are 


u= 2 Uy cos(n6) (2a) 


eo 
and V=>. v,s/n(n@&) , 
n= 


The necessary derivatives are 
~ 
f= — n*U, cos(n@) 


® 


00 (2b) 
Ve oe ra COSCNO) 
Substituting Equations 2 into Equation 1 we obtain 
Ep= EDM t EO Wu, +n Va} cos (ne), (3) 


The assumption that plane cross sections remain plane 


requires that the axial strain be 


€,= atb(ntr)cos(e), (4) 


where a and b are constants to be determined. 
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Because there are no external forces or moments applied 


to the tube So dA = 0 )) 
[yc dA=0 (6) 
[x@dA= 9 (7) 
where y=(nt+r)cos@ and dA=(nt+r)drdé@. (8) 


Equation 7 is automatically satisfied as a consequence 
of symmetry. 

With the assumption that radial stress is zero, the 
stress condition in the tube is that of plane stress. 


The equations of plane stress are 








Oz = ee [en te _- Ct VE (8) | (9) 
ol =, eerie Co “)E (a) | ; (10) 


where Eis the thermal strain. Expressed as a Fourier 


series E(e) =E,+> e cos( nem @ne 
f 
T oT” 
where €, = Af [andr do 
De 
and Ey = Lf [amdr cos(nO@)dO +; n=1,2,3,°*° 
Qo Tr 


The thermal strain is an even function as a consequence of 
the temperature being symmetric about the y-axis. In 


T() 
general E(a)=f a(t) dT - Using the above expression 
Tz 


late 





for axial stress and integrating Equations 5 and 6 permits 


determining a and b in terms of displacement: 











a= (1+yE,- LEMs (12) 
b = £ [c+v€, 7 (13) 
2 x a 
where § = I~ oon y Balt jaion, and §= /+ Fx 


The integration is performed in Appendix II. 
The elastic strain energy for the case of plane stress 


as discussed on page 4 is 


U = tf {lee EC), +[eg-E@] og} dA (14) 


per unit length of tube. 
Substituting Equations 9 and 10 for stress into Equation 


4, multiplying and collecting terms, we get 


ao = flee ref ravge,-2(- EGE, Ed +2(14E(8) [dA (25) 
4 


Integrating Equation 15, we get the following expression: 


. - oy LZ 2 z 
WV xs 20 +qSb> +4 [2E ur +2 Gun +20 glint, +0, )] 


= 2 atau. te b(uiry)] -2(i4y)[20&, +B % bE, | 


as alo) [2 BYE, +5. (au +1 Ve dE, | +2(weer+e 62], (16) 


See Appendix II for details of the integration. 
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Applying the theorem of least work as discussed pre- 
viously, we are now in a position to evaluate the coeffi- 
clients u, and v, by taking the partial derivative of X. | 
with respect to each coefficient and equating the result 
to zero. 


For — coefficient u, 


= 4a(£ ) + Lig + tal 4 a v£ ) 
AC VEE) - tCinyBE. 5 (17) 


Therefore, “=P r, a where a 3 sas : (18) 
Horeeceriicacncs u, and’ vy, , 


= 466 ++.lu,ty,)°+ aX Bb(u,+y) ~ 20 Br bE, 
ality) (u,+ ye, i 


where only those terms that contain u, or v, have been 
written. Note that b contains the sum ( u, +v, ). It is 


not possible to evaluate u, and v, separately because 


om equals OX . However, it is sufficient to obtain their 
‘ VY, 


sum because that is the only form in which u, and v, appear 


as can be seen from Equation 19. Taking 594— and equating 


Ut Y) 
to zero: 
ss b 2tu4 Vv, ){1- REDE: /- =O. 
Bo = tabby aah » ae udfi- 78'] -200& [1- B1)-0 (29) 
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Solving for (u,+v, ) we get 


_ Be 
(ut+vJ=2=¥e€, where p= Serie (21) 


Forvall other coefficients, u, and v, , n7l : 


> Ds u,tn Mn Or -(1+V)P, be E,= 0; (22) 
dx _ - 
and 57 = Pn Unt nV, -C+¥) 6 E,= O (23) 
yields nvy=(ItWGEE,- Pn &n. (24) 


a 
Substituting Equation 24 into 22 yields Cra JU,= 0 : 


a 
But CF,-2 )#O _ , therefore, 
Ly, = O 5 Vy, = GtV) Ip ba 2 (254a0n 


The radial and tangential displacements are now 
known as functions of the thermal strain, and thus tempera- 
ture. Substituting Equations 18 and 21 into Equations 12 


and 13 respectively, yields: 


a= ? E, where p = (+v-FY) (26) 
and b= p e, where fp = & [1+ ¥-4)], (27) 
iC g 


Thus the axial strain is 


E,= PE, t+ £)PE cose). (28) 


14 


| 
& & 





Substituting Equations 18, 21, and 25 into 3 gives the 


tangential strain, 


e,= YE, +WE,cos (8) +42 E, cos (n6) | (29) 


Remembering that E(a) = E, +E cos(@) +2 En cos (n@) 


then 2 E, cos (n 6) = E (6) ese ey cos O and, 
é€,=l(¥,-1)-vIE, +l(Y-)-v]E cos G +(1+v)E(8) . (30) 


Substituting these expressionsfor strain into the stress 


formulas, Equations 9 and 10, determines the stress: 





cen -£€(0)+ EMG) YG-) 16-v IE 1E-)+vd-Ye VIF fecose} (321 


{= 





and dg, = E {+(e +1G-) +E") + vp £IE,<os 6} (32) 


i-p* 


This concludes the stress analysis for the symmetric case. 


B. Anti-Symmetric Case. 
If the temperature distribution is anti-symmetric 
with respect to9=0, the thermal strain is anti-symmetric. 


Expressed as a Fourier series 


E(o)=>. E, sin(ne) , (33) 


A=/ 


Jee 


a ia 
_ 
— « @e ee 18 A 


AY eB a 


- '¢ » © Gees 


A. <q © kh 


= 


a _ @= 4 eae «& — —_— =: 


— te a 





T (6) 


where €,= [ [stodr siIn(nade. 


Note that there is no constant term in the Fourier series. 
It 1s seen that this is similar to the symmetric case 

rotated 90 degrees and we can apply some of the results of 

the previous case. 

The first harmonic is the only one contributing to axial 

strain, so the axis of rotation is the y-axis and the 


longitudinal strain can be written as 
6,=@2,+ 6,4 sinG@, (34) 


where a,is automatically zero because of anti-symmetry. 
The radial displacement will be an odd function and the 
tangential displacement an even function of the angular 


coordinate. Expressed as Fourier series they are 


“@) =D 4 : OD (35) 
and v0) = 2. Vn COSCO) . (36) 


Note that v contains no constant term. <A constant term 
would imply a rotation of the tube. 


Substituting Equations 35 and 26 into 3, we obtain 


eg = £> {lit Er* vi, - nV }sin (ne). (37) 
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Wetting —v, = Ay; yields 
eel = ST f,% oF 
6= £ > {lt bo 0]a, +d, fsin(na) . (38) 
o 7 


The expressions for € and €,are of the same form as 
in the symmetric case except that sin(n@) is substituted 
for cos(N@). The only integrals that do not vanish when 
integrated over the interval 0 to 27 are Isin? (nde. 
These integrals have the same values as Joos (n@)dé Or 
the interval. Thus by making the substitution Y=-¥4, and 
writing sin(m@) for cos(n@), the solution of the symmetric 


case is applicable and the results are as follows: 








é,= (+7 PE sin 6 (39) 
e,= Y-)-VJE sing + (i+ WE (@) (4.0) 
eo) at ~~ IE-) +vY-1) +0-v) +9 £ JE,sin 6 (41) 
o, = a lw-1) % V(P-tvd £ JE sin é, (42) 


The stress analysis is now complete for the anti-symmetric 


case. 
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C. General Case. 


For a general temperature distribution the thermal 
strain can be expressed as the sum of the symmetric and 


anti-symmetric cases, an even and an odd function. 
co CO 
Thus €(@)=€, +2 €, cos(n@) +2 E,sin(n@) > (43) 
a complete Fourier series. The general solution is as 
follows: €,= PE ++ 2 )¢ [E,cos@ +E siné] (44) 
Ep = 14-9 -VJE, HG -) -v]le,cosar+E smo] +G+vyEO@) — (45). 


= -EE@)+ FL f1G-1) +oG-+-vIIE, 
+I(G-) +vG%-+0-vy +B £J[E,cos 0 +E, sel} (46) 


O = & 1-0 + (B= NE, H1G-Y+WG-) +0 £ |lE,cosa- E, sine} }(47) 


4-yr 





The stress and strain resulting from a general temperature 


distribution is now known. 


D. Simplified Equations for Stress and Strain. 

It is instructive to consider the magnitude of — , 
a; g , and Y in order to simplify the equations of 
stress and strain. The stress and strain equations (Equa- 

tions 45, 46 and 47) contain the terms (f-/ ie (Dp -1 we 
(Y-/ ), and (Y-/ ) which vanish if the ¢'s and W's are 


1& 





unity. Table I shows that for tubes of practical wall 


thickness this is very nearly so. 


TABLE I 


The Effect of Wall Thickness on Stress and Strain 










foryv= 0.296 


substituting unity for P, ; v ; p, and Y, 
Equations 44, 45, 46, and 7 become: 


ee E. +(i+% Ecos + €, Sin Q) 
€,= C+WE(o) ~ v(E, +€)c0S 8 +E, sin @) 
Oz = -E{E(6) =[€.+01+ Ne, COS O+E, sin 9)}} 


— Z 
0; = Ge. (E, cos 9 +E,sin 8) | 


These equations should be sufficiently accurate for most 


design applications. 


i 


Ty, [veer | a.00o1e 00005 





(48) 


(49) 


(50) 


(51) 





4. Bowing of Tube. 
A tube subjected to the assumed temperature distribu- 
tion will bow into the are of a circle, or if restrained 
develop stresses equal to those required to straighten the 
bowed tube. This is of practical concern in pipe design 
and has been discussed by Flieder, Smith, and Wetmore [6]. 
From the analysis of the bending of an initially straight 
beam (7] it is known that the curvature, K ; is » where 
c is the distance from the neutral axis and€ is the strain 
resulting from bending. The axial strain on the mid-thick- 
ness curve of the tube is E,=2,E,+£ (Ecos +E, sin), obtained 
by taking r=O in Equation 44. The $ E, term denotes a 
general elongation of the tube and does not contribute to 
bowing. The d e cos O term results in bowing in the yze- 
plane. The strain, € , at 9=0 is b&,, at a distance r, 
from the neutral axis. Thus, Kye= oe is the curvature 
in the yz-plane. The PE, sm@ term results in bowing in 
the xz~plane. The strain at G=% igs PE, and as before the 
distance to the neutral axis is *% . Therefore, Kys* 26 


So 


is the curvature in the xz-plane. The curvature in general, 


K= EEE (52) 


lo 


is the vector sum of % and Ayz, and the plane of bowing 


olen 


is given by = tan ; (53) 
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measured clockwise from the y-axis. Again it is possible 
to set P, =1 without significant error so that 


ea (eae 


a : (54) 


ei: 





5. Results, Design Equations. 

The analysis leads to the following simplified equa- 
tions which can be used to design thin tubing systems 
subjected to steady temperature distributions that are in- 
dependent of tube length. 

1. Axial strain 


€,2E&,+(/+€ )(E,cos6 +E, s/n 6) (48) 
2. Tangential strain 


E,=(1+WE 6) - Y(E,t+E,cos 6 + E sing) (49) 


3. Axial stress 
op -E fe (0) -[é.+ (1+ aaa PIE, cosoré, sinel]} (50) 


4. Tangential stress 


og = PEE (Ecos @ +E, sino) (51) 
5. Bowing curvature 
: = “1 £ 
K= rte et Op=tan = (54) &(53) 
¥ T() “ 
ven E(a) =f «tr) dT =E. +> [e..cos@a)té, sin(n@)| 
Tease n=l 
aw TC) 
E, = tf f«(ndT dO 
© TB 
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6. Discussion of Results. 

When thermal strain is expanded in a Fourier series, 
the tangential stress and bowing as given by Equations 51 
and 54 are functions of the first harmonic only, and the 
axial stress as given by Equation 50 is a function of the 
second and higher harmonics. To interpret the physical 
meaning of this, consider the case where EQ)=E, cosat€,cosza, 
When the thermal strain is projected onto a yz-plane,the 
first harmonic transforms to a linear function and the 
second harmonic transforms to an even function. It is 
known [8 | that linear thermal strain does not cause stress 
in a free beam, but causes bowing, and that an even thermal 
strain function causes stress but not bowing. Thus, the 
results are those expected from superimposing elementary 
solutions for simple beams. 

Equations 44 thru 47 include a first approximation for 
the effect of wall thickness; however, the basic assumptions 
of this analysis, particularly that the temperature does 
not vary through the wall, suggest that the results be 
confined to thin walled tubes. Accordingly, Equations 48 
thru 51 and 54 are suggested as being the appropriate ones 
to use for actual design and analysis. These results 
afford a slight extension of those obtained by Flieder, 
Loria, and Smith in the following respects: (1) they give 


a variation of axial stress through the wall thickness and 
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(2) they give a value for tangential stress. 

The present results and those of Flieder, Loria, and 
Smith are comparable for the case where E, vw, and & are 
constant. When these conditions are imposed, Equation 50 
is identical at the wall mid-thickness to the result devel- 
oped by Flieder, Loria, and Smith. Note that 9, is 
independent of Poisson's ratio. Equation 50 gives a minor 
variation of oO through the wall; they did not consider 
any radial variation of 03 . Flieder, Loria, and Smith do 
not give a result for tangential stress. Equation 51 shows 
that tangential stress is zero at the tube wall mid-thickness 
and that it is small for tubes with slight bowing. Equation 
51 also shows that tangential stress vanishes if Poisson's 
ratio is zero. 

Comparison of the results for bowing show that there 
is a very slight effect due to finite wall thickness and 
non-zero Poisson's ratio as given by Equation 52. However, 
the simplified equation, Equation 54, which is here recom- 
mended for purposes of design and analysis, is essentially 
identical to that given by Flieder, Loria, and Smith. 

Essentially, we have shown that influences neglected 
in the analysis of Flieder, Loria, and Smith have only the 
very slightest effect on the final result. However, the 


expression of the thermal strain as a Fourier series has 
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Simplified the final form of the equations and aids in the 
visualization of the effect that a given temperature dis- 


tribution may be expected to have on a tube. 
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Simplified the final form of the equations and aids in the 
visualization of the effect that a given temperature dis- 


tribution may be expected to have on a tube. 
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APPENDIX I 
Derivation of Expression for Tangential Strain 

The temperature change from point to point around the 
circumference of the tube results in thermal strains. There 
is assumed to be no temperature gradient across the wall 
thickness of the tube, thus the thermal strain is indepen- 
dent of radius. Superimposed upon this strain is another 
strain which varies with radius and is caused by the change 
in curvature of the tube wall. 

To derive an expression for that portion of the strain 
independent of radius, consider two points A and P, close 
together, on the mid-thickness curve. When the temperature 
changes, the points will move to new locations A, and P, as 


shown in Figure 5. 





J 


Figure 5. Tangential Strain at Mid-thickness. 
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Lét £ and y be unit vectors in the i and j directions 


respectively. The bar denotes a vector. 


Then OA,= vi+(Qt+u)y , (1la-1) 
6Q =(rrur+du)[sin@O)z + cos (d9); | ; ieee 
QP =(vt+dv)lcosdQi + sin(d@jl , (1a-3) 
ana A, RP = OG + QP — OA, | (la-4) 


Because 09 is a small angle, sin(d9) =d9 and cos(d6) =1. 
Then, AP =Kn+Wde+ dud@+dvjz +[du-vde-dvd6l; . (Gees) 


Neglecting infinitesimals of higher order, 
AP =Un+u)dd+dv]z +[du+vd6] 7 (la-6) 
2. Zz 2 e ya ° 2 2 z 4 
and 1A,Pl=[(ntwd6* +2 (ctu)yd0"+a7d0-2evdesvd0l. (1-7) 


The strain, <9) »on the mid-thickness curve is defined 
as the change in length per unit initial length of the curve 
as the tube wall deflects. For small angles the arc length 


and the cord length can be taken as equal. Therefore, the 


strain is _ [A Pl-%d9 . 
Go/~ 4do. * (la-8) 


Substituting the value of IA,2l yields 
i 
[Cora)de+2lnu)ide+ cde —2avde"+ v*d6"]*— r dé? 


ae eee eee _ 
Eg nde (la-9) 


Again neglecting infinitesimals of higher order, the strain 
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La S + 
becomes . /+=(utv) ee 
Eg,/ . rl f Yr J = 2 (1la=-10) 
Oo 
az @ a e a za At v | 
Since pe Cos | , [1+2 C+) | A “Ke . (1a-11) 
Therefore, ° 
~ UstVv 


on the mid-thickness curve. 

The second component of strain is caused by the change 
in curvature of the tube wall. Because the tube is "thin" 
the theory for bending of straight beams can be used. The 


strain due to bending is 
a Ane (1a-13) 


| 


where r = distance from mid-thickness and 
AK 


Curve. 


ll 


change in curvature of tube wall mid-thickness 


To find the change in curvature consider a point P 
which after a change in temperature is at a new location 


E 


/ as shown in Figure 6. 
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Figure 6. Displacement of a Point on Mid-thickness 
Curve. 


The point is located by the vector R= (r+) E, + ve, (la-14) 
where @r and @€, are unit vectors in the radial and tan- 
gential directions respectively. Taking derivatives; 

R= 46+ (6+WE,+ VE, - ve, 
(4-v)é, +(K+UutvV) ey (la-15) 


and R= (4-V)E. +(A-V)&++We,-GKt+u+vée, 


=(k-2y-u-njJé, +(V+24a-V&, , (la-16) 
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The curvature after deflection is[{3] 


u-vila+2a—y)—-Ch -2V-Uu-KXr tury 


: la-l 
[G-v—* + (ctutv)*]% — 


K = 
Noting that r, greatly exceeds u, v and their derivatives, 


we get the approximation 


bes n+h@ut3y—-au) 
oe [i+ 2 (q+) ]% 


ee a) ee (La-18) 
“ee 


The change in curvature is the final curvature minus the 
initial curvature. The tube is initially round so the 


initial curvature is 1/y. Thus, 


Sa (la-19) 


| 
oa 5 f 


The strain due to bending is then 


=-Gt4)r (1-20) 


ADE 
9,2 


The tangential strain is the sum of the two components 


developed above: 





“tv kt Ur 
Ee, = Ky “7 As r (1) 
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Solving Equation 2a-4 for a, we have 


a= C4W&,- I~ £, ]¥¥ (2a-5) 
and @= ((+VJE,— — a 
where lea ! eer. (12) 


Equation 6. 
JG YdA=0 (6) 


y =(te+")cos @ and dA=(GtN drd@ (8) 


Substituting Equations 8 and 2a-l into 6 and integrating 


over the same wr as before: 


*O+ v2 + Sie $ E (n=! Jugtn V, bcosin8) cose 





~O+W2E, cos(n@)cos Aldo = O, (2a-6) 
Ww 5 me=yi 
Note that [cas m0) costno)d0 = 0; msn - (2a-7) 


Henceforth in the analysis, terms which vanish in this way 


will not be written. Continuing the integration, 





“ry (u4 Vix - (4V)E, (n+r) Idr=0 (2a-8) 
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Let $= I+ Fre and (eo = sae 
then b= A —[O)E, - vGew J. (13) 


Equation io. 


aT = fl len +€g +2VEz Eg -A +E) &y (15) 
L, 15 


-2(nvJE@e,+20+vE@) IdA 


The integral is the sum of the following terms: 
% nar 


I S@ dA = cradr fla b(t) cos* 6] dG 


% 
= J cae latwa‘+TbCn+r)*I]dr 
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where (= i+ pha ln) and Q,= I+ en (n*-1) 
am Save, é,dA= “eonil IE )ate+ (n+ntury,) beast] dé 
ve [Ceol E)adst(nerlury)bidr 
= BLT ty a goes: +(ut+y,) b6 
Iz | -2lev)E(@)e,dA= -2 (+) 4 Sadia + bG+JE cos*O ld 
=-R conf a§,+ b(n+E, en dr 


we: ee) faaE.+ br, BE, ] 


% 
ays [-26nE@) dA = 2C/+V)} ( coda E fh eii-luen a fE,cos’ ha dg 
A ; eae 
% 
= a MoE > Sir EG Nurn¥eddr 
2, 


= aes 2£u€, +> Quen VE, | 


a 
VI [eG»)@? de = 2 cry andr g E, cos *(ne) de 


=2ti+y)r5etl2 EP 4D ef | 
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Summing all terms: 


ACE = Xa2a% he's Lele Fur +> Guz +200, Mtn hl] 
+BY [eau F tury)bb]-2levea£,+ bb rE, | 


=o Cd) [eFUE, +2. (A u,tn VE y, | 


+2d+vl2eEs wy, S| (16) 
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TEMPERATURE, ‘°F 


APPENDIX III 
Example 

- The following example will illustrate the method of 
solution and order of magnitude of results. <A 302 stain- 
less steel pipe 10 inches in diameter, with a wall thick- 
ness of 0.5 inches is subjected to the temperature 
distribution shown in Figure 8. Assuming the pipe is free 
from restraint, what is the stress pattern and how much 


will the pipe bow? 
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Figure 8. Temperature of Pipe. 
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The first step is to determine the thermal strain for 
each temperature. The base temperature in this case is 
taken as 68°F because the data, see Figure 9, is given in 
that form. A harmonic analysis is then performed to deter- 
mine eon a , and E, » The calculations and results are 


given in Table II and paragraphs below. 





-400 -300 -200 -/00 O 100 200 300 400 $500 
TEMPERATURE, °F 


Figure 9. Thermal Expansion of AISI 302 Stainless Steel 


From T. F. Durham, R. M. McClintock, and R. P. Reed, 
"Cryogenic Materials Data Handbook", National 
Bureau of Standards, PB Report 171809, 1961. 
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TABLE II 
Harmonic Analysis 


9 T(Q) EG) a cos 
radians x /O? 


EC) sin 
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/ 49% 
ae) = - 128.3 X/o 
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E, = 7 5 El@)cos6 =/22.4x10- 


10 


~~ 


E, = cS a E@Msin@= 32.3xX oa 


The bowing of the pipe is 


_ ex Z = z _9¢ 
K= as ol = VGaed Sy 6220) oy Stax Ym 
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at Op= tan & = fan’ 323 = /%8 
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The tangential stress is given by 
E ~— 
09 = aE (E cos +E, sinG), 


therefore, on the mid-thickness curve, r=0, @=0. 
Letting E=30x10 psi and Y =0.296, 


we. GaecGoD . 27. 
I-bY 17-9087 P gots 


At the outside diameter r=t/2=0.25 inches. 


G 
05 = a 9.73> (122.4 COSO + 32.3 sin@O ). 
The results are given in Table III. 


The axial stress is given by 


Oz= -F{E(@) -[E.+ OF ao ME, cos + a Se 6)]} 
é | <> 
2 -30° $s@ -fni28 +i + py 224 Cos 9 432.3°sin 6)]} : 


fe 
On the mid-thickness curve r=O0 and [ / + Kr O-v*) | =], 
At the outside diameter r=0.25 and (1-4 Sap } =-1 coo 
Axial stress was calculated for each position and the 


results are given in Table III and also shown in Figure 10. 
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TABLE III 


Tangential and Axial Stress 


at r=0.25 at r=0 at r=0.25 






radians 
0.0 











AXIAL _Fe 0. 
STRESS 45 2:25. 





TANGENTIAL Lr=0.25 _ 
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Figure 10. Stress Pattern of Pipe. 
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The result of the first assumption, that 6, and Trg 
are negligible, is that equilibrium cannot be satisfied. 
The values of these stresses are now estimated using the 
conditions of the example. A section is cut out as shown 
in the free body diagram, Figure 11. These are the posi- 


tions of maximum and zero tangential stress. 





st 
Tre 


Figure ll. Shear Stress 
The stress distribution, og , causes a net counter-clock- 
wise moment of 25.7 in-lbs. To restore equilibriun, 

/ 
stresses 0, and @,, are required. Summing moments about 
the point P, 

Ti Te 
rO- OS(s) pao 

The tangential stress, 6, , is also equal to 10.3 psi. 

An indication of the magnitude of the radial stress 
required to maintain equilibrium can be obtained by inves- 
tigating the condition at 9@=195°, the cold side of the pipe. 
There the axial stress is about 10,000 psi and will be 


taken as constant. Figure 12 shows the action of axial and 
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radial stress. 





Figure 12. Radial Stress Component Caused by 
Axial Stress. 


The radius of curvature, R, was obtained from the 


bowing. Summing forces in the radial direction, 
Rome oC 7 Os ao O 


Therefore, 0-= Ope eee ) = 1.27 psi, the radial 


stress caused by axial stress. 
The maximum radial stress component resulting from 
tangential stress can be estimated by considering the stress 


condition at QO = OR as shown in Figure 13. 
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Figure 13. Radial Stress Caused by Tangential 
Stress. 


Taking the sum of the forces in the radial direction 
we get, 
2q.%~de-2 £(4q+G')dé 
es a 4G +6') 


- a25 [600 = ; 
O-= 225 ($8 + /0) 275 psc . 


The total radial stress is the sum of the component 
caused by axial stress and the component caused by tangen- 
tial stress ( O-= 427+ 72.75 =Q psd ). 

It is seen from the above calculations that 9d, and 
Cre are small when compared with other stresses, and there- 


fore, neglecting them is realistic. 
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